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A  Locally  Most  Poi-.'erful  Rank  Test 
For  the  Location  Parameter  Of  a  Double  Exponential  Distribution 

1.  SUTTMARY.   The  locally  most  po'-erful  rank  test  (L.M.P.R.T) 
for  the  location  parameter  of  the  double  exponential  c.d.f. 


L(x,0)  = 


|e(^-®)    X  < 


l-|e^^-)   x>0 


is  examined.   This  test  was  originally  found  by  Blrnbauni  [2]  whose 
unpublished  derivation  is  given  for  completeness  in  Section  2. 
A  remark  of  Ansari  and  Bradley  [1],  T,ihich  asserts  that  the  Chernoff- 
Savage  theorem  [5]  holds  under  certain  weaker  conditions  is  used 
to  prove  the  asymptotic  normality  of  Birnbaum's  test.   Coraparing 
this  test  with  the  likelihood  ratio  test  and  making  use  of 
Pittraan's  definition  of  asymptotic  relative  efficiency  (A.R.E.) 
we  find  that  Birnbgum's  test  is  asymptotically  efficient.   The 
A.R.E.  of  the  latter  to  the  likelihood  ratio  test  is  one  for 
syiTETietric  distribution  and  othsrijise  is  shoTin  to  vary 
between  zero  and  infinity.   The  latter  result  is  interesting  in 
contrast  with  the  situation  for  the  normal  distribution  for  which 
the  A.R.E.  of  the  L.M.P.R.T.  (the  Fisher-Yates)  to  the  likelihood 
ratio  test  (the  "t"  test)  is  always  greater  than  or  equal  to  one. 

2 .   Introduction  -  The  Birnbaum  Test 

Let  x^,X2,...,x^  ,  Y-^»72,"*>Y^     be  N  =  n^  +  n^  independent 
random  variables,  the  x.'s  having  unknown  absolutely  continuous 

c.d.f.   P^  and  the  y. 's  having  unknown  absolutely  continuous  c.d.f. 
X  1 

*  Alternatively,  a  theorem  of  Hajek  under  slightly  different 
assumptions  proves  the  same  result. 


t-ojij:. w  v:v;: J:\;    I^i  tS^'Li:^^.. ^M^'!?^. 


-I'.i;  rriij'' 


P  •  ^'^e  consider  the  problem  of  testing 


against 


Hq  :  P^  =  L(x,5^q),  Py  =  L(y,j2^Q) 


H^  :  P^  =  L(x,Q),  Py  =  L(y,0), 


where  -oo  <;^<jz^q<6<co   and  we  use  only  the  ranks 
Let 


'Nj 


0 


if  the  j   smallest  of  the  pooled  sample 
of  N  observations  is  an  x 

otherwise  • 


Capon  [t|.]  hfls  shovn  that,  under  regularity  conditions  ^'hich  are 
met  here,  the  L.M.P.R.T.  is  of  the  form 

0    if  B(z.^-,...,z_)  <  C 


(x^,..,x^,y^,...,y^) 


1    if  B(z^^,...,z^^^^)  >C 


\K    i^  ^^^jv-'^™^  =^ 


M' 


where 


%  =  Bj^(^)  =  B(z.^^,...,z^^^)  =  ^^  4^  Sn\i 


^Ni  =  \^^     !  #•  1°S  ^'  (^-'i") 


®  =  ^o' 


X  (x,Q)  =  L'(x,Q)  =  I  e"'^"®'    -O)  <  x  <  oo 

and  V,  is  the  i^^  smallest  of  the  N  observations.   The  notation 
E^  /   indicates  expectation  is  taken  under  the  assumption  that 
the  x«s  and  y's  are  distributed  according  to  -C   (x,e)  and  -t!(y,©) 
respectively. 


['fe  find 
E 


^0^0  V"^  '''    ■''''^'' 


^  -  ^fird 


where 


1 

i.  1 


:J     <     ^^ 


?hus 


4u  =  P  U  or  "^ore  x^  >  j!^q] 


N-i  or  less  x,  <  j^  '^ 
i    0' 


=  i:(  1 )( 


l!^  (  ^  )(  |)^ 


)  -  1 


Letting  B-j  =  —  ^  ^'\r*^Mi»  ^  st^-tistic  originally  described  in  [2], 

we  see  that  the  Birnbaura  test,  a  linear  combination  of  binomial 
coefficients,  is  the  L.M.P.R.T.  for  the  location  parameter  of  the 
exponential • 


3.   Asymptotic  Normality  of  B^ 

Capon  [l\.]   has  obtained  from  Theorem  1  of  Chernoff  and 
.Savage  [$]    and  a  simple  eztension  of  their  Theorem  2,  that  a 
L.M.P.R.T.  has,  asymptotically,  a   normal  distribution  if  a  certain 
smoothness  condition  is  satisfied.   In  our  case  the  condition  is 


J(^)(H)|  = 


d^'-)j^(H) 


dH 


(i) 


<  K(H(1-H))"^"  2  "^  ^ 


for  i  =  0,1,2   and  for  some  5  >  0  where  k  is  a  constant  and 


'^  ■  -^^M' 


where  -i 

^  ('-1   H  <  •! 

I  1    H  >  i 

This  is  cleorly  sstisfied  for  all  H  except  for  the  point  H  =  -^ 

where  the  derivative  does  not  exist.   However,  Ansari  and  Bradley 
[1]  have  remarked  that  the  conclusion  of  the  Chernoff -Savage 
theorem  holds  even  with  such  a  discontinuity.   They  assert  that 
this  follows  immediately  upon  inspection  of  the  proof,  (the  present 
author  has  not  checked  this)i'"  Prom  this  we  conclude  that  for 


absolutely  continuous  G_{x) 


1     x2 


lim  Prob^  .'  .I—2A iL  <  t  I   =   /    (2Tt)  ^     B     ^     dx 

L  ^e^^  N''      ^     -OD 


00 

'^o^^V   =  I  'i}   -I  ^V^^^^  +^a^(x)]dGg(x) 


=  1  -  2G^(x") 

.V,  n^        n^  -^ 

and  where  x   is  that  x  for  'hich  |r-  G^Cx)  +  jf   '^q(^)  =  p 


Also 


I  J   JG^(x)(l-G^(y))j/^(Hg^(x)J^(Hg^(y))dGg(x)dGQ(y) 


00  <  X  <  y  <  00 

I 

0.(x)(l-G.(y))J,  (H  ,(x)J,  (H.  ,(y))dG,(x)dO,(y)  '(  , 


12         f  f       - -•■    .  .-• - 1 


^1 

-CO  <x<y<oo 


In  a  recent  article,       Jaroslav  Hajek  has  pointed  out 
that  asymptotic  normality  can  be  proved  under  other  regularity 
conditions  '-'hich  we  moy  nx!A   here  without  subst-intialljr  changing 
any  of  the  results.  (Dec.  1962  Annals  of  Math.  Stat.) 


where 

Of  course  the  conclusion  is  not  valid  if  :tv  -(B  )  =  0.   In 

2        2 
particular,  under  the  null  hypothesis  '-^  /  ,/  =  ^^   is  given  by 

n,N    p         }      ^  1  2 


In  the  non-null  case  we  calculate 


2    ^^    ^    _  ^n^        I      [    f      G^(x)(l-G^(y))gg(x)gg(y) 


^^S/^n)  =  T    i  ^'  '/ ."  "7     "'  "    dJ#Kx)H<^y)) 


.oo<x<y<co  Hg^(x)H^^(y) 


np         r   ^    G^(x)(l-G^(y))g/(x)g.(y)  1 

""  TT      J  J  ^ ^ dJL(H(x))dJj.(H(y))| 


1       ■'    ■■      ,,  H      (x)H      (y) 

-co  <x<y<oo  '^  ^ 


2,    ■--. 


n^     (  G^(x')(l-G^(x"))g^(x")        n^     Gg(x")(l-G^(x"))gJ(x   )| 
(H,.(x"))2  '  (H,.(x,))2 


Of  course  this  calculation  restricts  G  to  only  those  c.d.f.'s 
such  that  Hg,(x')  is  not  zero  where  ^  <  x"  <  ©  • 

h •   Asymptot ic  Relative  Efficiency 
i''e  now  use  the  Pittman  [8][7]  criterion  for  the  A.R.E.  of  two 
sequences  of  tests  N  =  |  ^ ^  ''^  and  ^>'"  =  |  w "  ]     .   Consider  the 

sequence  of  alternatives  /  .t  =  '^tvt  -  Kt  =  ^^N"  2  where  k  is  a  non- 
0,j  -  0      n,  n 

zero  constant  -^ ~  = =•  and  /^  =  Q  -  0.  Let  lim  — -  = 

e  -Q      ^2      — '  2 

n-T  ... 

=  lim  — T-  =  r.   Suppose  ^-^  and  w  have  the  same  size.   The  A.R.E. 

^2 

of  ^f  to  ^'"  is  defined  to  be   ,,^^  —   where  N"  is  the  sample 

size  of  the  second  test  required  to  achieve  the  same  poi.jer  for  a 


.    r 


given  alternative  iJiich  ".,^  achieves  with  respect  to  the  same 
alternatives  when  using  a  sample  of  M  observations* 

Consider  such  a  sequence  of  alternatives  x^iith  /\^   =  0  -  / 
and  a  sequence  of  statistics  ^-^  =  )  ^^ n  ^'  ^^'^   ^'''"  =  '•■   ^''-t  '^^    satisfying 
the  following  three  conditions  in  the  neighborhood  of  /^  =  0(^=0= 

t      1  .2 


(a)    JXp,fV%!V    ]-.   i     .-^^   ax 

n^  n„\2 


exists  and  is  independent  of  k. 

E  (P)  is  called  by  Pittrnan  the  efficacy  of  ^'^   at  P  and  P 
refers  to  the  distribution  under  which  expectations  are  taken. 

Pittrnan  has  sho^^'n  that  if  (a),  (b)  and  (c)  are  satisfied  for 

H  =  ^  v^  \     and  V*  =  I  T.r^  j  then 

"K-   E  (  ^  ) 

N->oo   N   g  ..(pj     W,W  ^^' 

(where  E,^^  v''^^^^    ^^   *^®  A.R.E.    OP  i-'  to  w'") ,    if  E,,,-;KP)    t^  0   and 

lim       ^1        _,.      "'l,^ 

N->oo      ^     -  lim  -.  -  r   . 
2  n 


To  find  E„  (G)  we  need 


-k%^^V 


/x 


/\=0 


7\  ]  ^  -  2G^(x  )  ; 


A=o 


2n^   >. 


N    &tJ-  li 


G  (x) 


\i=f 


where  x  is  the  median  of  g/  (x).   To  see  this  recall 


=  V->  --1^  ^V''' 


for  some  /?$,  ©  >  0  > 
Thus 


2n. 


2H^^(x)  -2G^(x)  =--^/^^G^(x)|^^. 


2n, 


1  -  2G, 


(x""')  =  -  -r/  A  -r-G  (x'"") 


[l=fi 


Differentiation  yields 


A.=o 


2n„   ^     ^ 


as  was  indicated,   ^'e  have  then 


V°'=w_^V^'l^o'' 


Ghernoff  and  Savage  have  pointed  out  that  the  asjrmptotic 


t'...-:'  ■' 


8 

efficiency  of  a  sequence  of  tests  can  be  established  by  means 

of  comparison  with  a  likelihood  ratio  test.   They  have  also  shown 
asymptotically  that  for  C  ,  compared  to  the  t  test,  C„  is  ali-mys 
more  efficient  except  at  the  normal  distribution"  where  the  A.R.E. 
is  unity.   Capon  has  shown  that  under  sviitable  regularity  condi- 
tions, not  satisfied  by  the  double  exponential  distribution,  the 
A.P..S.  of  the  L.M.P.R.T.  to  the  likelihood  ratio  test  is  unity. 
A  natural  question  one  may  ask  is  whether  the  Chernof f-Savage 
result  holds  for  all  other  distributions.   A  negative  answer  is 
obtained  by  investigating  tlie  double  exponential  case.   That  is, 
we  shall  show  that  the  A.^.S.  of  Birnbaum' s  test  to  the  likelihood 
ratio  test  L   of  ti-e  double  e::ponential  ranges  from  0  to  oo  and  is 
unity  for  L.   ''e  find 

L„  =  —   j  number  of  x.    >  ^^  -   number  of  x.  <  ^q] 
-   — —   [  number  of  y,  >  jz^q  -  number  of  y .  <  j6q^ 

2     '  -        2 

=  1 -[■  nmnber  of  x.  >  ^^  i   -  1  +  — 

n,    I  1   ^0  •        n2 


j  number  of  y.  ^  S^q] 


=  2tS'^0>  -  ^'^o'l 


where  L^{»)    is  the  sample  curaulative  distribution  function  of  a 
sample  of  n  random  observations  on  f;.   To  calculate  the  efficiency 
of  L  at  G  we  find 


Chemoff  and  Savage  considered  distributions  satisfying  certain 
regularity  conditions  and  Hodges  and  Lehmann  [6]  extended  to 
arbitrary  distributions. 


-i    '*     &"* 


■f'.ii    ciJ-^X 


^Q/h^""^    =  f   -I9  1°S  l'(x,. 


(gg(x)  -  g^(x))dx 


=    iO-^)      I      ^    log  l'(x,Q)I      ^r-Sn^^M   ~  d^ 
, ,-^ 

for  some  ^,  9  <  ^  <  ^  . 

Thus  assuming  g  satisfies  regularity  conditions  sufficient  to 

allow  interchf^nge  of  differentiation  and  integration  yields 


^Wh^""^ 


A=o 


^ 


[X  °p. 


(x)  dx 


-.  s'^'L- 


~G    (x) 


+  —  G    (x) 


2c^ 


S 


^ 


G.,(j2^n)  +  t::  G.,(-oo  )  .     +--—  G,,(co) 


M^  ^ti^^O 


<^[i    \i 


[i    [i 


Also 


(r   |g)  =  (  ^  +  ^)    . 


^0     N 


A  particularly  simple  case  results  when  we  consider  translation 
or  scale  alternatives,  i.e.  Gg^Cx)  of  the  form  G(x-9)  or  G(|), 
Then  combining  our  results  yields 


Clearly  when  G  is  the  double  exponential  the  A.R.E.  is  unity 


as  it  is  whenever  /^_,  and  the  median  of  g  coincide.   However  g  may 


r:::-f 


10 


be  chosen  so  thqt  E   n  (G)  ranges  betveen  zero  and  infinity. 


For  example,  if 


go(^)  = 


2^' 


Qx/27i 


then  the  median  x  =  p.   Thus,  at  i^S-  =  1, 


E.   R  (G) 


6e 


\   Ig^(P) 


and  as  P  gets  sufficiently  close  to  0,  the  A.R.E.  becomes  arbitrar- 
ily small. 

Similarly,  let 

'       "? 

X  >  0 


~2(a+l)   ® 


z^y-) 


(ae-Q  x)   gQx 
2(a+l) 


X  <  0 


where  a  >  0  and  6  >  0.   Then  the  median  x  =  0  and  at  ^q  =  1 


( a+1 ) e 


-1 


As  a  gets  sufficiently  close  to  0,  the  A.R.E.  becomes  arbitrarily 
large. 

5 .   A  remrk  on  the  r6l?,tionship  botveen  3.j  and  Mood's  M  test 
It  is  interesting  to  compare  the  Birnbaum  test  and  Mood's 

test  M.-  [3]   For  the  latter;  asymptotic  normality  has  been 

established  [3]  with  mean  n-G  (x")  and  variance 


n^n^F^ix^^Ml-F^^x^O) 


io•^    Z''y- 


11 


Calculation  similar  to  those  for  B,.  yield 


E^(G) 


^  C-  (x)!   , 


Thus  M  and  B  have  A.R.E.  unity  for  all  distributions  satisfying 
the  above  regularity  conditions.   In  particular,  this  gives  an 
asymptotic  justification  for  use  of  Mood's  test:  I'hen  the  double 
exponential  shf^pe  is  true.  Mood's  test  has  A.R.E.  unity  (B-^  seems 
preferable  under  this  assumption  since  it  has  the  exact  locally 
best  property  ?s  well  as  A.R.E.  unity). 
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